Heat conduction and Wiedemann-Pranz Law in disordered Luttinger Liquids 
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We consider heat transport in a Luttinger liquid (LL) with weak disorder and study the Lorenz 
number for this system. We start at a high-T regime, and calculate both the electrical and thermal 
conductivities using a memory function approach. The resulting Lorenz number L is independent 
of T but depends explicitly on the LL exponents. Lowering T, however, allows for a renormalization 
of the LL exponents from their bare values by disorder, causing a violation of the Wiedemann-Franz 
law. Finally, we extend the discussion to quantum wire systems and study the wire size dependence 
of the Lorenz number. 



PACS Numbers: 71.10.Pm, 72.15.Eb, 72.15.Nj 

The Wiedemann-Franz (WF) law states that for elec- 
trons elastically scattered by impurities the ratio of the 
thermal (k) and electrical (cc) conductivities normalized 
by T is a universal number Lq — n/ocT — 7r^/3e^ called 
the Lorenz number. This law is obeyed for electrons in 
conventional metals [Q , and inclusion of electronic inter- 
actions that lead to the standard Landau Fermi-liquid 
(FL) behavior does not modify Lq Hi, as their effect on 
K and (Tc scales away when the particle-particle scattering 
rate is much smaller than the impurity scattering rate. 

A fundamental question arises whether the WF law 
survives in strongly correlated electron systems with non- 
FL ground states. Interest in this question has been re- 
cently revived by an observation of violation of WF law 
in high- Tc cuprates |^ which are widely believed to be 
non-FL systems. An important and tractable example of 
a non-FL is a one-dimensional (Id) Luttinger liquid (LL) 
with separate charge (p) and spin (cr) bosonic collective 
excitations instead of fermionic quasiparticles. Heat 
being carried by both spin and charge excitations, there 
is a potential for violation of WF law in a LL. However, a 
potential complication in the LL case is that even a weak 
disorder can cause Anderson localization (7 - 11 1 . 

Experimentally, thermal transport in (quasi-) Id has 
been investigated in spin-chain, spin-Peierls, and 
spin-ladder materials There are also some early 

measurements of heat conduction in quasi-ld organic 
conductors Theoretically, although crdT) in the 

LLs has been extensively discussed [p^-p^ with mem- 
ory function methods |l9[ |, the thermal conductivity and 
the WF law have only been addressed in two contexts. 
First, Kane and Fisher (KF) examined spinless LLs 
with a single strong impurity for repulsive electronic in- 
teractions, and found a T-independent but nonuniversal 
Lorenz number L^^^^ . The generalization of this problem 
to the bulk disorder case is however not straightforward 
pl| (as discussed below). Second, Fazio et al. and 
Krive p3| have studied thermal transport in a mesoscopic 
Id wire connected to leads. They showed that even for 
a clean wire the contact with leads causes violation of 
the WF law except at sufficiently low T. However, the 
effect of backward scattering by finite-density impurities 
has not been considered in the latter case. 



In this letter, we first examine the WF law for an 
infinite LL with a nonzero concentration of impurities. 
Unlike KF, we calculate k{T) and adT) in the high-T 
regime where impurity scattering can be treated by the 
memory function method pg| , p^ . The resulting Lorenz 
number is independent of T but depends on the LL ex- 
ponents. It differs from the one obtained by KF. As T 
goes down, we show that the Lorenz number acquires T 
dependence as the consequence of instability towards An- 
derson localized state. Such a violation of the WF law is 
in sharp contrast to the higher-d FL situation. Finally 
we examine the effects of contact in a realistic quantum 
wire measurement and discuss the wire-size dependence 
of the Lorenz number. 

Conductivities from Memory function method. A LL 
of spin-1/2 fermions in the presence of impurities is de- 
scribed by the following Hamiltonian 



H = Hll+H,,^p = j dxHix), (1) 

i—p,a 

Himp = J dxri{x)dx(l)p 

+ [ dx{^e'^^^^'='>cos[V2c^„{x)]+ H.c.l, (3) 



where Ki are the LL exponents, Ui the velocities of the ex- 
citations, a the lattice constant, and Iii{x) and (f)i{x) are 
canonical momenta and coordinates, respectively, ^(cc) 
in Eq. (^) is the 2kp component of the impurity po- 
tential, assumed to have a Gaussian distribution with 
zero average, ^{x)^*{x') — D^S{x — x'). rj{x) is the for- 
ward scattering component. The charge and heat cur- 
rents Jc = J dxj^c{x) and Jq = J dxj'g^x) are obtained 
from the continuity equations dxJ^c{x) + dtn{x) = and 
dxJqix) +dt'H{x) — 0, where n{x) = —^dx(t)p{x) is the 
electron density operator. We find 

Jc{x) = V2upKpnp{x), (4) 
Jq{x) = ^ulllpd,(j)p - uln^d^cj), + JP{x), (5) 
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where Jf^ix) = ^/2upKpr]{x)Ylp{x). We note that the 
transformation (j)p{x) — (f>p{x) — ^P^Mf- v{^')dx' ehmi- 
nates rj from both i?imp and Jq , implying a null effect of 
r](x) on (Tc and k. Thus we set ri(x) = hearafter. Eqs. 
(||) and (||) also show that the heat and charge current 
operators are quadratic and linear in boson operators, 
respectively. This is crucial in resulting in deviation of 
the Lorenz number from the universal number at high T 
even for spinless fermion case as shown below. 

In the absence of impurities, both Jc and Jq are 
conserved currents. In fact, Jq is (up to a prefactor) 
the total momentum of the system, and its conserva- 
tion results from translational invariance. With impu- 
rities, neither Jc and Jq are conserved and finite con- 
ductivities can be expected At sufficiently high 
T, quantum effects are cut off by inelastic thermal pro- 
cesses, and disorder scattering can be treated within the 
Born approximation. This yields conductivities to lead- 
ing order in D^^. A convenient formalism to do this 
is the memory function method PqHTtUT^ , in which the 
finite-frequency electrical and thermal conductivities are 
expressed in terms of the memory functions Mj(uj) = 
Xj' [m;F,))^ - {{F,;F,))^^o] {j = c, Q) 



as 



iiiLo)^tXQT-'[Lo + MQ{uj)]-\ 



(6) 
(7) 



with Xj = T J^^'^ dX{Jj{Oy Jj{iX)) the current static 
susceptibilities, Fj{x) = [J'j{x),H], and ({Fj;Fj))^ — 

- J dxdx' dt e*'^* {[Fj{x,t),Fj{x' ,0)]) h- To leading 
order in Dj, both Xc and XQ take their pure LL values. 



7 UpKp, XQ ^ f (up + u„) T\ and 
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U 



2ttT 



A-Kt] T 



(8) 



where Kt = Kp + K^^ B{x,y) is the Euler Beta function, 
ac = l-Kt, UQ ^ ac-2, Ac{uj) = i^pUp/vr, and Aq{uj) = 

tt[{K^/4: + {cj/2ttTY] [Kpu^ + K„ul)/[Kt{Kt + 1)]. In- 
serting Eqs. (|i) into Eqs. (^[|7|), and after some straight- 
forward algebra, we find that the dc conductivities be- 
come 



2-Kt 



k{T) 



2e^T{Kt) a (2TTaT\ 
7rT^Kt/2)^ \ ) 
27r^Tl + Kt jup + u^)^ r{Kt) 
9 Kt ulKp + ulK,V^{Kt/2 



= e^lAT), (9) 



f2TTaT\ 



2-Kt 



oca(r) /ei(r), 



(10) 



where V = {2D / TTu^){ua / u p)^ '' is a dimensionless dis- 
order parameter pT| ], r(a;) the Gamma function, Cy{T) cx 
T the specific heat of a pure LL, and lc\{T) the elastic 
mean free path Both ac and n exhibit power-law 
behavior in T with nonuniversal exponents. 



Lorenz number at high T. From Eqs. (^ and (|T^), 
one has 



L 



(up + u^f 1 + Kt 



cjcT 9e2 u^Kp + ulKa Kt 



{ID 



Kp, Kcr, Up and are bare intrinsic parameters, so L is 
independent of T, and a generalized WF law is obeyed. 
In the noninteracting case, Kp = Kg- ~ I, Up = Ug- = vp, 
the universal Lorenz number L = Lq is recovered from 
Eq. (|ll|). 

The result ([ll] ) can be reduced to the spinless LL case, 
by making Up — Ua- = u, Kp = K^ — K. We obtain 



L' = (7rV9e2)(X-2^ 2/^-1). 



(12) 



It follows that L' < Lq for an attractive interaction case 
K > 1, as a result of tendency towards a supercon- 
ducting state with high electrical conduction but poor 
thermal conduction; While for a repulsion case K < 1, 
L' > Lq, indicating tendency towards a weakly pinned 
charge density wave (CDW) state with better heat trans- 
port than charge transport. It is instructive to compare 
L' with L(kf) ^ (^2/g2)(^2 _^ 2K)-^ obtained by KF 
for a single strong impurity at K < 1 We see that 
L' - L(kf) ^ 2{K - iy/9K^{K + 2) > 0. From physical 
point of view, L(^^) and L' can be roughly understood 
as the results for a strongly pinned and a weakly pinned 
CDW phase, respectively, so that we should expect a 
larger L^^^'^ than L' . This puzzle can be resolved by 
noting that some Hamiltonian terms causing heat but no 
charge conduction arc neglected in the derivation of the 
tunneling Hamiltonian used to calculate L^^^^ ||2^. In- 
cluding these terms explicitly for K = 1/2 indeed leads 
to |§ L(^^) = |L(^^) > L'. 

We would also like to remark that the single strong 
impurity case considered by KF corresponds to a low 
energy fixed point which cannot be reached from the 
limit of weak impurity scatterers at high concentra- 
tion c we consider |gl|. The KF fixed point is real- 
ized when the renormalized impurity strength tsil) = 



tsiO) becomes of the order of the high energy 



3(2-2^)/ 

cutoff W ~ u/a. This requires the renormalized length 
ae'-' — a[W/tB{0)]^^^'^^'^^^ to remain still much smaller 
than the inter-impurity distance 1/c. In the high impu- 
rity concentration limit, c — *■ 00, is — *■ with D — ct^ 
fixed, obviously such regime cannot be observed. 

Lorenz number with decreasing T. The enhanced 
quantum interference effects, which are responsible for 
occurrence of the Anderson localization at T < Tioc — 
Ua/hoc {hoc = aV^^^^"^^^*^ 0)i lead to renormaliza- 
tion of Kp, K„ and T). Now we study the influence of 
such renormalization on L for T still much larger than 
Tioc, so that a perturbative renormalization group (RG) 
theory can be used. We first neglect renormalization of 
Kp and Kcr by V. The RG flow equation determining the 
renor malization of V reads [0 dV{l)/dl = (3 - Kt)V{l), 
which leads to 
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V{l)IV = e^^-^'^K (13) 

It is clear that Kt — ?> defines a metal-insulator transition 
(MIT) line: for Kt < 3, !?(/) grows exponentially with 
renormalized length, and the Anderson localization will 
take place eventually; Whereas for Kt > 3, !?(/) flows to 
zero, corresponding to a delocalized phase. By dcflning 
a T-dependent scale r(r) from ae''(^) = u„/{2nT)_^ 
hh{T) {hh the thermal length), we see that Eq. ( |l^ ) 
yields {a/nV){2naT/u„)'^-^' = ae^'<^'^y[nV{l*)], which 
allows for a simple RG interpretation for CTc and k in Eqs. 
(p[]l0|): They can be obtained by applying the RG flow 
up to the scale 1*{T) at which T becomes order of the 
energy cutoff TiUc/a. At this scale, thermal effects sup- 
press quantum effects, and the use of the memory func- 
tion method is justified. One can then perform a memory 
function calculation in which the bare parameters T) and 
a are replaced by !?(/*) and ae' , respectively. 

A similar analysis can be made for inclusion of renor- 
malization of Kp and by V. We start from the bare 
parameters Kp{{)), Ka{Q) and £'^(0) and follow their RG 
flow, being described by Eqs. (3.4) in Ref. |ll[], up to the 
scale 1*{T). Since D^{1*) < 1 stiU holds for T > Tioc, 
but the quantum phase coherence has started to lose, 
Eqs. (p[]lO|) are valid again |2j], but with the bare Kp 
and Kcr replaced by the renormalized ones, Kp{l*) and 
Ka-{1*), respectively. The resulting RG- improved L ac- 
quires T dependence through Kp and Kcr, and the gen- 
eralized WF law holding at high T breaks down. This is 
in contrast to a higher-d FL case, in which such a renor- 
malization is always negligible . 
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FIG. 1. Lorenz number as a function of temperature 
close to the MIT line for spinless fermions. W ^ u/a 
is a high energy cutoff. The initial conditions are 
K{0) = 3/2, D(0) = 10"^ 

We proceed with a discussion of the spinless fermion 
case for simplicity. The above analysis suggests that L' in 
Eq. (|l^ as a function of the initial parameter K becomes 

i(T) ^ (^V9e^) {[Air)]-2 + 2[Kin]-'}, (14) 

where K{1*) can be obtained from the RG flow equation 



0, dK{l)/dl = -K'^{l)V(l)/2. For a generic initial K 
away from the MIT line, the RG flow line in the K-V 
plane is almost vertical, and the deviation of L from its 
high-T value is in order of T) which is insignificant unless 
T becomes close to Tioc. Alternatively, considering the 
second order terms in D in tJc and n one also obtains a 
T-dependent power-law correction term to L{T) with the 
exponent predicted by the RG method. 

However, in the vicinity of the MIT line, the RG flow 
is no longer vertical, leading to a faster variation of L 
with T. \K{1*) — K\ is found to be proportional to 
- yJV{l*) - V, indicating that L{T) shown in Eq. (|l|) is 
systematic. L(T) for this case is shown in Fig. ^ We em- 
phasize that such a behavior is not expected in a higher- 
d FL case, in which the Lorenz number is constant even 
close to the MIT point |]. 

Lorenz number in dirty quantum wires. An experi- 
mental realization for the measurement of the transport 
properties of a LL is through a flnite size wire connected 
to leads or a carbon nanotube ||2^ . The theory pre- 
sented here for a bulk material needs to be modified, as 
the boundary condition imposed by the contact is shown 
P?! , ^ to drastically influence the transport through the 
wire. Charge crosses the contact in terms of pieces of 
fractional charges |^^, whereas heat is carried by plas- 
mon modes and they cross the contact like waves passing 
through barriers [p2| . The Lorenz number will obtain a 
nontrivial wire size dependence as shown below. For sim- 
plicity we only discuss about spinless LL wires, and leave 
the discussion of the spinful case for a future publication 
pgf . Besides, we follow Ref. and model the leads as 
Id noninteracting FLs. 

We first briefiy recall the clean wire case. The elec- 
trical conductance G is purely from the contact and is 
quantized [|^, G — Go = e^/27r (for h ~ 1), which im- 
plies a perfect transmission of the total charge after a 
long time (dc) measurement The mismatch of the 

plasmon wave velocities at the two sides of the contact 
leads to T- dep endent thermal conductance /Cth normal- 
ized by T [g2[, which, at low temperatures T <^u/d id 
the wire size), asymptotically approaches to a constant 
/Cth/T ~ 7r/6 H. The WF law is typically violated for 
a generic T and is restored at T ^ u/d. 

Let us now turn to the case of a wire with impurities. 
We consider the case of hoc being much larger than d 
and/or /th so that we can neglect Anderson localization. 
The total electrical resistance is the sum of the contact 
resistance G^^ and the wire resistance G~l^^, resulting in 
the following total conductance, 

G = {G^'+G^J-\ (15) 

Gwiro depends on the relative values of the two length 
scales d and Zth(T). 

1) At d ^ Ith (T ^ u/d), a simple Ohm law holds 
for the electrical conductance, i.e., Gwiro = '^{T)/d = 
Go2TTL\/d with cr(T) and /d being found in Eq. (|^). 
Eq. ( |l5|) reduces, for d ^ ^oi(T), to Maslov's perturbation 
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result |g], G ~ Go[l - d/(27r/oi)], and, for d » lei{T), 
crosses to the Drude formula. In this regime, the ther- 
mal conductance of even a clean wire, /Cth, from Eq. (9) 
of Ref. 11^, depends on the details of the barrier, so 
that the thermal conductance of the dirty wire becomes 
a rather non-universal function of disorder, wire size, and 
temperature [ p9| . 

2) At d < /th (T < u/d), Eq. (9) of Ref. || indicates 
that heat is carried by plasmon modes of wavelength 
much larger than the distance between the contacts, and 
thermal conductivity of the clean wire becomes universal. 
In this regime, we can make use of the T 4^ u/d equiv- 
alence |3^,Q to obtain an explicit d-dependent Lorenz 
number. According to Ogata and Fukuyama, the finite 
size effect can be taken into account within the memory 
function formalism by replacing to in Eqs. (^) and by 
i2u/d. This leads to 



G(d) 
/Cth(d) = 



1 



TT 2u + 2uCciu,K)f3{d)' 
■KuT 1 



3 2u + 2uCt^{u,K)(3{dy 



(16) 
(17) 



with I3{d) = {d/hocf-^ , and Cc(m, KVa^d Cc{u, K) be- 
ing easily found by comparing Eqs. (|6|) and (^^ with 
Eqs. (^) and (|l^). In writing down Eq. ( p^ we have 
carefully taken into account the screening of the electric 
field coming from the leads js^. We recover a universal 
Lorenz number in this regime in the absence of impuri- 
ties in agreement with [ p2[ . In the presence of impurities, 
Eqs. dlO) and (|l^) immediately yield 



L{d) ~ io 



l+Cc{u,K)l3{d) 
l+Cth{u,K)(3{d)' 



(18) 



which is a function of the wire size since Cc{u,K) ^ 
Cth{u, K) in the presence of interactions. Thus, devia- 
tion from WF law is obtained in finite size systems with 
impurities even for T <C u/d. 

Conclusion. We have investigated the WF law in a dis- 
ordered LL system. At high T where thermal effects cut 
off Anderson localization, the Lorenz number is constant 
in T, and, in the spinless fermion case, smaller (larger) 
than the universal Lq for attraction (repulsion) between 
the fermions. Its dependence on the LL exponents is dif- 
ferent from the one obtained in a single strong impurity 
case [Q. As T goes down, Anderson localization effects 
induce an interplay between disorder and electronic in- 
teractions, which is responsible for the violation of the 
WF law. When a dirty Id wire is connected to leads we 
find the Lorenz number becomes a function of the wire 
length. 
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